CONSTRUCTIVE PERTURBATION THEORY FOR MATRICES 
WITH DEGENERATE EIGENVALUES* 
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Abstract. Let A(e) be an analytic square matrix and Aq an eigenvalue of A{0) of multiplicity 
m > 1. Then under the generic condition, J|: det (A/ — A (e)) | ;^)_(q^jj) 7^ 0, we prove that 
the Jordan normal form of A (0) corresponding to the eigenvalue Aq consists of a single m X m 
Jordan block, the perturbed eigenvalues near Aq and their eigenvectors can be represented by a 
single convergent Puiseux series containing only powers of s^/"*, and there are explicit recursive 
formulas to compute all the Puiseux series coefficients from just the derivatives of A (e) at the origin. 
Using these recursive formulas we calculate the series coefficients up to the second order and list 
them for quick reference. This paper gives, under a generic condition, explicit recursive formulas 
to compute the perturbed eigenvalues and eigenvectors for non-sclfadjoint analytic perturbations of 
matrices with degenerate eigenvalues. 
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1. Introduction. Consider an analytic square matrix A (e) and its unperturbed 
matrix A (0) with a degenerate eigenvalue Aq- A fundamental problem in the analytic 
perturbation theory of non-self adjoint matrices is the determination of the perturbed 
eigenvalues near Aq along with their corresponding eigenvectors of the matrix A (e) 
near e — 0. More specifically, let A (e) be a matrix-valued function having a range in 
£"nxn^ the set of 71 X n matrices with complex entries, such that its matrix elements are 
analytic functions of e in a neighborhood of the origin. Let Aq be an eigenvalue of the 
matrix A (0) with algebraic multiplicity m > 1 (which may be degenerate, i.e., m > 1). 
Then in this situation, it is well known [U §6.1.7], [21 §11.1.2] that for sufficiently small 
e all the perturbed eigenvalues near Aq, called the Ag-group, and their corresponding 
eigenvectors may be represented as a collection of convergent Puiseux series, i.e., 
Taylor series in a fractional power of e. What is not well known, however, is how 
we compute these Puiseux series when A (e) is a non-selfadjoint analytic perturbation 
and Aq is a degenerate eigenvalue of ^ (0). There are sources on the subject like [1], 
[5], [5], [TJ §7.4], and [3 §32], but it was found that there lacked explicit formulas, 
recursive or otherwise, to compute the series coefficients beyond the first order terms. 
Thus, the fundamental problem that this paper addresses is actually two-fold. First, 
find a method to determine how many Puiseux series there are that represent the 
Ao-group and their eigenvectors along with the fractional power of e that is associated 
with each. And second, find explicit recursive formulas to compute all the series 
coefRcients. 

Main Results. In this paper we will show that under the generic condition 

^det(A/-A(£))|(^,^^(„,^^^0, (1.1) 

the fundamental problems mentioned above can be solved. In particular, we prove 
Theorem 12.11 and Theorem 13.11 which together state that when this condition is true 
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then the Jordan normal form of A (0) corresponding to the eigenvalue Aq consists of 
a single m x m Jordan block and for e sufficiently small there is exactly one Puiseux 
series for the Ao-group and one for their corresponding eigenvectors both of which 
contain only powers of e™. Theorem 13. II further goes on to state that the Ao-group 
and their corresponding eigenvectors can be represented by a single convergent series 
of the form 

A, (e)=Ao + ^A, (C^^) 

oc ^ 

Xh (e) = xo + 'Y^ Xk (c'^e™) 
fe=i 

for h = 0, m — 1 where ^ = Ai ^ 0, and xo is an eigenvector of ^ (0) 

corresponding to the eigenvalue Aq. More importantly though, this theorem gives 
explicit recursive formulas that allows us to determine the A/c's and x/c's from just the 
derivatives of A (e) at e = 0. Using these recursive formulas, we compute the leading 
Puiseux series coefhcients up to the second order and list them in Corollary [331 

The results of this paper depend essentially on the assumption that the generic 
condition is true. In theory or practice, however, it may not be a straightforward 
task to show its true. In order to deal with this issue we provide Theorem l2.1l This 
theorem gives alternative statements that are equivalent to the generic condition (|l.ip . 

The key to all of our results is the study of the characteristic equation for the 
analytic matrix A(e) under the generic condition (II. By an application of the 
implicit function theorem, we are able to derive the functional relation between the 
eigenvalues and the perturbation parameter. This leads to the implication that the 
Jordan normal form of the unperturbed matrix A (0) corresponding to the eigenvalue 
Aq is a single m x m Jordan block. From this, we are able to use the method 
of undetermined coefhcients to get explicit recursive formulas for determining the 
Puiseux series coefficients. 

We want to take a moment here to show how the results of this paper can be used 
to determine the Puiseux series coefficients up to the second order for the degenerate 
case m > 2. To begin we may assume, because of Theorem 12.11 that the Jordan 
normal form of the unperturbed matrix A (0) corresponding to the eigenvalue Aq 
consists of a single m x m Jordan block. This implies we may put A (0) into the 
Jordan normal form [S] §6.5: The Jordan Theorem] 



U-^A{0)U = 



Jm (Ao) 






Wo 



where (see notations at end of §1) Jm (Ao) is an m x to Jordan block corresponding 
to the eigenvalue Ao and Wq is the Jordan normal form for the rest of the spectrum. 
We next define the vectors ui,. . . , Um, as the first to columns of the matrix U, 

Ui :— Uci, 1 < i < m (1-3) 

(These vectors have the properties that ui is an eigenvector of A (0) corresponding 
to the eigenvalue Ao, they form a Jordan chain with generator Um, and are a basis 
for the algebraic eigenspace of A{0) corresponding to the eigenvalue Ao). Next we 
partition the matrix U~^A (0) U conformally to the blocks Jm (Ao) and Wq of the 
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matrix U ^A{0)U a.s such 



U^^A {0)U 
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(1.4) 



Then it follows from Theorem 12. II and Theorem 13. II that 

d_ 



f det(A/-^(e))|(e,A)=(o,Ao) 



/ SI 



- dct{XI-A{e))\(^c,x-) = (o,\o-) 



(1.5) 



And thus the generic condition is true if and only if am,i 7^ 0. This gives us an 
alternative method to determine whether the generic condition p.ip is true or not. 

Lets now assume that am.i 7^ and hence that the generic condition is true. 
Define / (e, A) := det {XI — A (e)). Then by Theorem 13.11 and Corollary 13.31 there is 
exactly one convergent Puiseux series for the perturbed eigenvalues near Ao and for 
their corresponding eigenvectors given by 



2 

A, (e) = Ao + Ai (C^^) + A2 (C^" ) + E (C'e- 

fe=3 

Xh (e) = xo + xi (c'^e^^ + xi [c^^e'^^ + ^ a;/c (c''e " 



fe=3 



(1.6) 
(1.7) 



for /i = 0, m — 1, where ^ 
order may be given by 



\ l/m 



A2 = 



\ m! d\ 

^rn— 1.1 ^ 

\rn — 2 



f (0,Ao) 



Furthermore, the series coefficients up to second 
\ 1/™ 



{K^'^'i^ (0, Ao) + Ai^ (0, Ao)) 



toA 



m-l f J_ 
1 \m\ 



1^(0, Ao! 



Xo = Ui, Xi = \iU2, 

-AA' (0) ui + X2U2, if TO = 2 

X2 = < ^ ^2 

A2U2 + AjU3, it TO > 2 



(1.8) 

(1.9) 

(1.10) 
(1.11) 



for any choice of the TOth root of a^.i and where A is given in 

The explicit recursive formulas for computing higher order terms, Afc and Xk, are 
given in Theorem 13. II using (|3.12p and (j3.13|l . 

Example. The following example may help to give a better idea of these results. 
Consider 



A is) 



i 
1 1 





' 2 





1 — ^ 


+ e 


2 





-1 




1 





_ 



(1.12) 
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Here Aq = is an eigenvalue of A (0) of algebraic multiplicity m — 2. We put 
A (0) into the Jordan normal form 
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-1 
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U-^A{0)U = 











,u = 
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, u-' = 


-1 
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1/2 J 
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1 
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-1 



so that Wo — 1/2. We next define the vectors ui, U2, as the first two columns of the 
matrix U, 





'1 




'l 


Ml := 







1 




1 




1 



Next we partition the matrix U ^ A (0) J7 conformally to the blocks J™ (Ao) and Wq 
of the matrix U^^A (0) U as such 



u-^a'{o)u = 
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-1 


" 




" 2 





-1 " 




"l 
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l" 




" * 




-1 
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-1 
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1 1 
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* * 





Here 02,1 — 1, ai_i — 0, and 02,2 — 1- Then 

^det(A/-A(£))|(,,^)^(o,Ao) 



1 — 02,: 



^det(AJ-A(£))[(,,;,)^(0,;,„) 



implying that the generic condition p.ip is true. Define / (e, A) := det (A/ — A (e)) = 
A'^ — 2A^e — ^A^ + Ae^ — ^Ae + + |e. Then there is exactly one convergent 
Puiseux series for the perturbed eigenvalues near Aq = and for their corresponding 
eigenvectors given by 



A,, (e) = Ao + Ai ((-!)'' e^) + A2 ((-1)" e^) + J] A, ((-1)" e^' 
(s) =xo + x, ((-if si) + X2 ({~lt +Y.XU ((-if e--\ 

k=3 

for ft, = 0, 1. Furthermore, the series coefficients up to second order may be given by 

Al = 1 = 



/a2,i 



1 _ ai,i + a2,2 

2 ^ 2 



f (0>Ao) 



^0, 



-(A;jTft(0,Ao) + Ai^(0,Ao)) 
Ai (^1^(0, Ao)) 





"1" 




T 


Xo = 







1 




1 




1 



X2 = -AA (0) Ml + A2M2 
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by choosing the positive square root of 02,1 — 1 and where A is given in 



Here 



A = U 


[ Jm (0)* 


















{Wo- 


- Ao/n-m) 


-i 














"1 


1 1" 




' 







1 


-1 


" 




" 3 


-1 


-2 







1 1 




1 







-1 


1 


1 




3 


-1 


-2 




1 


1 




. 


(1/2)"^ J 
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so that 



X2 — — AA (0) Ml + A2M2 
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Now compare this to the actual perturbed eigenvalues of our example (|1.12p near 
Aq = and their corresponding eigenvectors 



A,, (e) 



1 







-1)'' 


c- 2 




"1" 




"f 


(e) = 
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"f 
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fc=3 



\h (e) 



/c=3 



for /i = 0, 1. We see that indeed our formulas for the Puiseux series coefficients are 
correct up to the second order. 

Comparison to Known Results. The example above illustrates the broader 
goal of this paper, to help bridge the gap between the theory and computability in 
determining the perturbed eigenvalues and eigenvectors for analytic perturbations of 
non-self adjoint matrices. The most relevant papers [4], [5], [6] and books [U §7.4], 
[71 §32] do not bridge this gap completely. For instance, the paper [6] addresses 
linear perturbations and gives theorems ^64 Theorems 2.1 and 2.2] more general then 
ours under which one may compute the leading exponents and leading coefficients 
for both the eigenvalue and eigenvector perturbations. But unfortunately, like the 
other references, when it comes to computing higher order terms these sources either 
do not address the issue or do not address it thoroughly enough. For example, of 
all the works mentioned only the work of Trenogin and Vainberg [71 §32 pp. 413- 
418] addresses the issue of higher order terms for analytic perturbations. But their 
work on the determination of the higher order terms does not lead to explicit recursive 
formulas like those we give in Theorem l3.1l ?3 of this paper. Furthermore, the formula 
for the second order term [T] §32 p. 417] and the implicit recursive formulas for the 
higher order terms [7, §32 p. 418] of the perturbed eigenvalues contain misprints. 
In our paper. Theorem 13.11 gives explicit recursive formulas for all the higher order 
terms and Corollary [33] gives correct formulas for the leading terms up to the second 
order. These results are achieved by following a method more similar to Vishik and 
Lyusternik [!4, Appendix I]. 
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Overview. Section 2 deals witli tlie generic condition p.ip . We derive condi- 
tions that are equivalent the generic condition in Theorem 12.11 In §3 we give the 
main results of this paper in Theorem [23 on the determination of the Puiseux series 
with the recursive formulas for calculating the series coefficients. As a corollary we 
give explicitly the leading Puiseux series coefficients up to the second order. Section 
4 contains the proofs of the results in §2 and §3. 

Notation. Let C"^" be the set of all n x n matrices with complex entries and 
C"'*^ the set of all n x 1 column vectors with complex entries. For a G C, A G C"^", 
and X = \ai,i]"^i & C"^^ we denote by a*, A* , and x* , the complex conjugate of a, the 



conjugate transpose of A, and the 1 x n row vector x* := [ 
y € C"^^ we let {x,y) := x*y be the standard inner product. The matrix / e 



1.1 



n.l 



For X, 



is the identity matrix and its jth column is G 



The matrix /„ 



is the 



{n — to) X (n - 
A to be 



identity matrix. Define an to x to Jordan block with eigenvalue 



Jm(A) 



A 1 



When the matrix A (e) e C"^" is analytic near e = wc define Ak := -gA^'^'> (0). Let 

2. The Generic Condition. The following theorem, which is proved in §4, 
gives conditions which are equivalent to the generic one (jl.ip . 

Theorem 2.1. Let A{e) be a matrix-valued function having a range in C"^" 
such that its matrix elements are analytic functions of e in a neighborhood of the 
origin. Let Aq be an eigenvalue of the unperturbed matrix A (0) and denote by m its 
algebraic multiplicity. Then the following statements are equivalent: 

(i) The characteristic polynomial det (A/ — A{e)) has a simple zero with respect 
to e at X = Xo and e = 0, i.e., 

^det(A/-^(£))|(^,^^(„,^^)^0. 

(ii) The characteristic equation, det(A/ — A(e)) — 0, has a unique solution, 
e(A), in a neighborhood of X = Xq with e(Ao) = 0. This solution is an analytic 
function with a zero of order to ai A = Ao, i.e., 



d'^e (A) 



dA" 



d'"-^e (A) 



A=Ao 



dX' 



A=Ao 



0, 



d™£ (A) 



dx'- 



A=Ao 



^0. 



(iii) There exists a convergent Puiseux series 



Xh (e) = Ao + AiC^™ + E Afe (c'^e^Y , h 

k=2 ^ ' 



where C, — e™*, which gives all the solutions of the characteristic equation, for suf- 
ficiently small e and X sufficiently near Ao- Furthermore, the first order term is 



nonzero, i.e.. 



Ai ^0. 
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(iv) The Jordan normal form of A (0) corresponding to the eigenvalue Xq consists 
of a single mxm Jordan block and there exists an eigenvector Uq of A (0) corresponding 
to the eigenvalue Aq and an eigenvector vq of A (0)* corresponding to the eigenvalue 



Aq such that 



(vo,A' (0)uo) ^0. 



Remark 2.2. Under the same hypotheses, the statement that (XqI — A{e)) 
has a simple pole at e — and the eigenvalue Ao has geometric multiplicity one, is 
also equivalent to those above in Theorem \2.1l Unfortunately, the proof of this lies 
outside the scope of this paper so we can only hint at why it is true. Basically, the 
statement given at the beginning of this remark or those given in Theorem \2.1\ are 
both equivalent to the local Smith form ^ p. 394] of the matrix XqI ~ A{e) at e = 
being 

XqI ~A{e)=E (e) diag{e, 1, 1)F (e) 

where E{e), F {e) G C"X" gj-g analytic and invertible matrices in a neighborhood of 
the origin. 

3. Determination of the Puiseux Series and the Recursive Formulas 
for Calculating the Series. This section contains the main results of this paper 
presented below in Theorem 13. II To begin we give some preliminaries that needed to 
set up the theorem. Suppose that A (e) is a matrix- valued function having a range 
in C"^" with matrix elements that are analytic functions of e in a neighborhood of 
the origin and Aq is an eigenvalue of the unperturbed matrix A (0) with algebraic 
multiplicity m. Assume that the generic condition 



^ det (XI 

OS 



^(^))|(s.AMO,Ao)^0' 



is true. 

Now, by these assumptions, we may appeal to Theorem 12.11 (iv) and conclude 
that the Jordan canonical form of A{0) has only one mxm Jordan block associated 
with Aq. Hence there exists a invertible matrix U g C"^" such that 



U-'^A{0)U 



Jm (Ao) 






Wo \ 



(3.1) 



where Wq is a (n — m) x (n — m) matrix such that Aq is not one of its eigenvalues [31 
§6.5: The Jordan Theorem]. 

We define the vectors wi,. . . , Um, wi,. . . , G ^nxi ^s the first m columns of the 
matrix U and (C/^^) , respectively, i.e., 

(3.2) 



K:=U 



u, : = Uei, 1 < i 


< TO, 


: = 




1 < i < m. 








[ Jm (0)* 








(Wo- 


' Xolji-m) 



(3.3) 
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where {Wo — Ao/n-m) ^ exists since Aq is not an eigenvalue of Wo(for rnore details 
on the vectors p.2p and the matrix A see Appendix A). 

Next, we introduce the polynomials pj^i = pj,i (Ai, . . . , Xj-i+i) in Ai,. . . , Aj_i+i, 
for J > « > 0, as the expressions 

POM ■■ = 1, Pj,o 0, for j > 0, (3.4) 

^ i 

Pj,, (Ai, . . . , Aj_,+i) : = 22 U ^s.Jor j>i>0, 

siH hSi=J 

l<Se<3-t+l 



and the polynomials r^^ij^ = r^+i^ki^i: in Ai,..., A/, for fc, Z S N, as the 
expressions 



rfc+i,fc(Ai) :=0, rfc+i,fe(Ai, A,) ^ J] ^s,, for / > 1 (3.5) 

siH l-Sfc = fe+i '^"^ 

i<s5<; 

(for more details on these polynomials see Appendix B). 

With these preliminaries we can now state the main results of this paper. Proofs 
of these results are contained in the next section. 

Theorem 3.1. Let A{e) be a matrix-valued function having a range in C"^" 
such that its matrix elements are analytic functions of e in a neighborhood of the 
origin. Let Aq be an eigenvalue of the unperturbed matrix A (0) and denote by m its 
algebraic multiplicity. Suppose that the generic condition 

^det(A/-^(£))|(^^,^^(„^,^)^0, (3.6) 

is true. Then there is exactly one convergent Puiseux series for the Xg-group and for 
their corresponding eigenvectors which may be given by 



k=l 

for h = 0, m — 1, where C = e^' with 

f , . f det(A/-A(£))|(,,^)^(o,^„) 
\ = {Vrn,AiUi) = - ^^0- 

V J 

Furthermore, we can choose 

Xi = {v^,A,u, f"' , (3.9) 

for any fixed mth root of {vm, ^iWi) and the eigenvectors to satisfy the normalization 
conditions 

{vi,Xh{e)) = 1, h^O,...,m-l. (3.10) 



A,(e) = Ao + 5]A,(c\^) (3.7) 

fe=i 

Xh{e)^xo + f2^kU''e^y (3.8) 
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Consequently, under these conditions Xi, A2,. • • and Xq, Xi, ... are uniquely deter- 
mined and are given by the recursive formulas 



( 



l/r. 



1/rn 



f det(A/- A(£))|(,,;,)^(0,A o) 

^ dct(\I-A{s))\^,_x) = {o,Xo) 



(3.11) 



s-1 s-1 

-r„^ + s-l.,l^+T, E Pj.i I I'm-!-: 
i = j = i 



+ Z; Z; Pj,i V^^i, J2 



for 2 < s < m 
X (3.12) 



k^m + s-l-j-kr. 



for s > m 



'^Ps,iUi+i, ifO<s<m~l 



i=0 

s-m j L~sr"J 



(3.13) 



PsAUi+1 ~ J2 Pj,kA''+'^AiXs^j-im, if s > m 

i=0 j=0 k=0 1=1 



where [J denotes the floor function. 

Corollary 3.2. The calculation of the kth order terms, Xk and Xk, requires 
only the matrices Ag, . . . , A^ m+k-i ^ . 

Corollary 3.3. The coefficients of those Puiseux series up to second order are 
given by 



Ai 



{VjnjAiUj 



sl/m 



- A 



i(O.Ao) \ 
^.§^iO,Xo)) 

(^r+^(7^|g^(O.Ao) + Aijj^(O,Ao) + 10(O,Ao)) 

-(^r + ^(^|g^(O.Ao)+Ai^(0,Ao)) 
'"^r"'(iTl^(0:Ao)) 

{vi,{A2- AiAAi)ui) , ifm^l 

i/m > 1 ' 



, ifm = l 



if m > 1 



(Dm- 1 ,AiUi ) + («,„ ,Ai«2) 


lA. 



Xo = Ui, 



Xi 



X2 = 



— A^iiti, if m — 1 
X1U2, if m > I ' 

-AA2 + (AAif - AiA^^i^ ui, if m=l 
—KAiUi + X2U2, if ni — 2 
X2U2 + A1U3, if m,>2 



where / (e, A) := det (A/ - ^ (e)). 

4. Proofs. This section contains the proofs of the results of this paper. We 
begin by proving Theorem 1 2. II of §2 on conditions equivalent to the generic condition. 
We next follow this up with the proof of the main result of this paper Theorem 13.11 
We finish by proving its corollaries 13. 21 and [ 
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4.1. Proof of Theorem 12.11 Proof. To prove this theorem we wih prove the 
following chain of statements (i)^(ii)=>(iii)=^(iv)=^>(i). 

1. (i)=>(ii). Suppose (i) is true. Define / (e, A) := det {XI — A {e)) . 
Then / is an analytic function of (e, A) near (0, Aq) since the matrix elements of 
A (s) are analytic functions of e in a neighborhood of the origin and the determinant 
of a matrix is a polynomial in its matrix elements. Also we have / (0, Ao) = 
and ^ (0, Ao) ^ 0. Hence by the holomorphic implicit function theorem [9l §1.4 
Theorem 1.4.11] there exists a unique solution, £(A), in a neighborhood of A = Aq 
with e (Ao) = to the equation / (e, A) = 0, which is analytic at A = Ao. We now 
show that e (A) has a zero there of order m at A = Aq. For (e, A) near (0, Aq) we can 
expand / into a power series / (e, A) = J2i+j=o i jet>i ^ij^^ (-^ ^ Ao)"*. By hypothesis 
Aq is an eigenvalue of A (0) of algebraic multiplicity m hence A = Ao is a zero of 
order m for the characteristic equation of A (0). This means that d^f\dX^ (0, Ao) = 
for < i < m - 1 but a™/\9A" (0, Ao) ^ 0. Combining this with the fact that 
/ (0, Ao) - and f (0, Aq) ^ we get 



OO 



/(e,A) =aioe + aom(A-Ao)"'+ }^ a^.e^iX-XoY (4.1) 

i+j=2, ijefi 
(i,i)^{(Oj):j<m} 

where aio = §j- (0, Ao) ^ and aom = ;j^fv^ (0; Ao) ^ 0. Now by definition e (A) is 
an analytic function that satisfies / (e (A) , A) = in a neighborhood of A = Aq with 
£ (Ao) = 0. Let q be the order of the zero of e (A) at A = Aq so that for |A — Ao| << 1, 

£ (A) = £q (A — Ao)'' + O ^(A — Ao)^^^^ with £g 7^ and q > 1. Then using the series 

expansion of / (£, A) about (0, Ao) and e (A) about A = Ao together with the identity 
/ {e (A) , A) = for |A — Ao| << 1, we can solve for q and Sq finding that q = m and 



1 3"* dct(AJ-A(e)) 
m! OA'" 



(A,£)^(Ao,0) 2) 



f det(A/-A(£))|(^_^^^^^^_^/ 

Therefore we conclude that e (A) has a zero of order m at A = Aq, as desired. 

2. (ii)=>(iii). Suppose (ii) is true. The first part of proving (iii) involves 
invert e (A) near £ = and A — Xq. To do this we expand e (A) in a power series 
about A = Ao and find that e (A) — g{X)"^ where 



5(A) = (A - Ao) £m + ^ £fc(A-Ao 



k—m+1 



and we are taking any fixed branch of the rnth root that is analytic at £„. Notice 
that, for A in a small enough neighborhood of Ao, g is an analytic function, g (Ao) = 0, 
and (Ao) = Em™ =/= 0. This implies, by the inverse function theorem for analytic 
functions, that for A in a small enough neighborhood of Ao the analytic function g (A) 
has an analytic inverse g^^ (e) in a neighborhood of £ = with g^^ (0) — Xq. Define 

a multivalued function A (£), for sufficiently small £, by A (£) :— g^^ (^'" ) ^^^'"'^ 

£ m we mean all branches of the mth root of e. We know that g^^ is analytic at £ = 
so that for sufficiently small s the multivalued function A (£) is a Puiseux series. And 
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since %^(0) 



we have an expansion 



k 

A(e)=ff"^(£") =Ao + Ai£™ +^Aft(£™) . 

fe=2 

Now suppose that for fixed A sufficiently near Aq and for sufficiently small e that 
dot (XI — A (e)) = 0. We want to show that this implies A = A (e) for one of the 
branches of the mth root. Well we know by hypothesis we must have e = s{X). 
But as we know this implies that e = e (A) = 5(A)'" hence for some branch of the 
mth root, b„i{), we have bmis) = bm{g(Xy^) = g{X). But A is near enough to Aq 
and e is sufficiently small that we may apply the g^^ to both sides yielding A = 
g~^ (5(A)) = (6m(e)) — X (e), as desired. Furthermore, all the m branches A/, (e), 
h = 0, m — 1 of A (e) are given by taking all branches of the mth root of e so that 



,(e) = Ao + AiC'e^+^Afe (c^^) 



fe=2 



where ^ = e ^ ' and 



dg^ 
dX 



(Ao) 



(4.3) 



3. (iii)^(iv). Suppose (iii) is true. Define the function y (e) := Xq (ff™). 
Then y is analytic at e = and || (0) = Ai ^ 0. Also we have for e sufficiently small 
det {y (e) I — A (e™)) = 0. Consider the inverse of y (e), y^^ (A). It satisfies = 
det (y (y-i (A)) / - A{[y-^ (A)]")) = det (A/ - A ([y-i (A)]")) with y-' (Aq) = 0, 

(Ao) = X^\ Define g{X) := [y-^ (A)]". Then g has a zero of order m at Aq 
and det (XI — A (g (A))) = for A in a neighborhood of Aq. 

Now we consider the analytic matrix A {g (A)) — A/ in a neighborhood of A = Aq 
with the constant eigenvalue 0. Because is an analytic eigenvalue of it then there 
exists an analytic eigenvector, x (A), of A {g (A)) — XI corresponding to the eigenvalue 
in a neighborhood of Aq such that x (Aq) ^ 0. In particular, x (Aq) is an eigenvector 
of A (0) corresponding to the eigenvalue Aq. Thus by using the series expansions of 
g{X), A [g (A)) — XI, and x (A) about A = Aq we get 

Q={A{g{X))-XI)x{X) 
= (a (Ar" (A - Ao)'" + O ((A - Ao)'"+')) - (A - Aq) / - Aq/) x (A) 

= ({A (0) - Ao/) - / (A - Ao) + Ar"^ (0) (A - Ao)") {x (Ao) + 

g (Ao) (A - Ao) + • • • + ^ (Ao) (A - Ao)"-^ + ^ (Ao) (A - Ao)^ 



+0 ((A-Ao)"'+') 



dx 



{A (0) - Ao/) X (Ao) + [{A (0) - Ao/) — (Ao) - x (Ao) ) (A - Ao) + • 
+ ((A(0)-Ao/)^(Ao)-^ 



(Ao)) (A-Ao)"-' 
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+ [{A (0) - Ao/) ^ (Ao) - (Ao) + K"-— (0) X (Ao) j (A - Ao)^ 

+0 ((A-Aor+'). 
This implies that 



(A(0)-Ao/)x(Ao) = 0, 
(A(O)-Ao/) — (Ao) = a;(Ao), 



(^(0)-AoJ)^(Ao) = ^(Ao), 



(7™T rl"^~^T HA 

{A (0) - Ao/) ^ (Ao) = (Ao) - XT""— (0) x (Ao) . 



The first m equations imply that ^ _i (Aq) is a Jordan vector of degree m 



-i-i-i 



with (A (0) - Ao/)' 1^ (Ao) = (Ao), i = 0,l,...,m- 1 so that .t(Ao), 

(Ao), . . ., (Ao) is a Jordan chain of length m generated by j^^-T (Ao). 

Since the algebraic multiplicity of Ao for A (0) is m this implies that the there is 
a single m x m Jordan block corresponding to the eigenvalue Ao where wc can take 

^ ('^o) ' 31 (•^0) ' dX"^-^ (^0) ^ Jordan basis. Furthermore, it is a fact from linear 
algebra that the Jordan normal form of a matrix and that of its adjoint are the same 
except the eigenvalues arc replaced with their complex conjugates. Thus using this 
fact we can imply that A (0)* has a single m, x m Jordan block corresponding to the 
eigenvalue A^. Let vq be a eigenvalue of A (0)* corresponding to the eigenvalue Aq. 
Then 

^0, {A (0) - Ao/) ^ (Ao) j = i^{A (0) - Ao/)* ^0, ^ (Ao) 

0,^(Ao)j=0. 
Now by the {m + l)th equation then we have 

0=Uo,(A(0)-Ao/)™^(Ao) 



Vo, 



T(Ao)-Ar"^(0)x(Ao))) 



vo, (Ao)) - Ar" (vo, ^ (0) X (Ao)) 



We will show that wc can choose a ijq an eigenvalue of A (0)* corresponding to the 
eigenvalue Aq such that (^vq, (^o)^ =/= 0. To do this we extend x (Ao),^ (-^o), 

. . . , §^ (Ao) to a basis x (Ao), f (Ao), . . . , £^ (Ao), Vm+i, Vn, for C"xi. We 
form a n X n matrix U with this iDasis as a column vectors as such 



U 



d"^~^x 

X (Ao) I • • • \ ^yn-l (-^o) \ym+l\ • • • 12/r, 
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Then since x (Ao), ^ (Ao), 



' rfA"*-^ (Ao), Um+i, Un is a basis it forms a linear 
independent set and so U is invertible. Define V as tlie matrix 



v = {u-^) 

Then we have the biorthogonal relation 

V*U = I. 



We will show that the mth column of V is the vector we seek, i.e., vq — Vcm. = 
(?7^^) Cm is an eigenvector of A{0)* corresponding to the eigenvalue Ag such that 

(yo, S^T^ (Ao)) 7^ 0. First we find that (vq, (Ao)) ^ since 

-(Ao)] - [{U-^) e,n,Ue„,j 



Now we need only show that vq is an eigenvalue of A (0)* corresponding to the eigen- 
value An. To do this we first note that 



U-'^A{0)U ■ 



where Wo is a (n — m) x (n — m) matrix such that Ao is not one of its eigenvalues. 
Hence by Proposition lA.ll fii) it follows that vq = {U~^) em is an eigenvalue of ^ (0)* 
corresponding to the eigenvalue Ag. Hence 



Jm (Ao) 






Wo \ 







Wo, 



dA" 



dA 



-(Ao))-Ar"(wo,^(0)a;(Ao) 



dA 



= 1-Ar™ (^wo,^(0)x(Ao) 

implying that (wg, ^ (0) x (Ao)) = A™ ^ 0. Therefore we have shown that the 
Jordan normal form of A (0) corresponding to the eigenvalue Ao consists of a single 
m X m Jordan block and there exists an eigenvector u of A (0) corresponding to the 
eigenvalue Ao and an eigenvector w of A (0)* corresponding to the eigenvalue Aq such 



that 



(iv)=^>(i). Suppose (iv) is true. We begin by noting that since 
det (Ag/ - A (e)) = (-1)" det ((A (0) - Ag/) +A' {0)e + o (e)) 
= (-1)" det (^{A (0) - Ao/) + A (0) e) + o (e) . 



it suffices to show that 



de 



det 



((A(0)-Ag/)+A' iO)e) l^.^O. 



An expression for the derivative of det (ji^A (0) — Ag/) + A (0) e) can be found in 

[TD| Theorem 2.16]. We will prove from this expression that there exists a nonzero 
constant c such that 



de 



■ det 



((A (0) - Ao/) + A (0) e) I^^Q = c («o, A' (0) ^^o) ^ 
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which will complete the proof. 

Let A (0) — XqI = VSX* be a singular-value decomposition of the matrix A (0) — 
XqI where X, Y are unitary matrices and S = diag((Ti, . . . , f7„_i, cr„) with ai > 
. . . > (7,1-1 > cTn > 0(see [3i §5.7 Theorem 2]). Now since the Jordan normal form of 
A (0) corresponding to the eigenvalue Aq consists of a single Jordan block this implies 
that rank of A (0) — XqI is n — 1. This implies that ai > . . . > Cn-i > cr„ = 0, 
u = Xcn is an eigenvalue of A (0) corresponding to the eigenvalue Ao, and v = Ycn 
is an eigenvalue of A (0) corresponding to the eigenvalue Ag. The latter implies that 
since the Jordan normal form of A{0) {A{0)*) corresponding to the eigenvalue Aq 
(Aq) consists of a single Jordan block then there exists nonzero constants ci, C2 such 
that u = CiUq and v = C2Vq. 

Now using the result of [TUI Theorem 2.16] we find that 

d 



— det ({A (0) - Ao/) + A (0) e) = 



where 



Sn-i=det{YX*) cT^,■■■a,^_,dei((Y*A{{))Xy \. 

Here ( y*A (0) X\ is the matrix obtained from Y* A (0) X by removing rows 



■il ...4„_1 

and columns «i . . . in-i- But since ct„ = and 



{y*A{Q)x\ =elY*A{Q) 

V / l...(n — 1) 

[v,X {Q)u 



C2C1 {vq,A (0)wo) 7^ 







then 



Sn-i^det{YX*) o-,, •••CTj„_, det ( (0)X 

l<il<---<i„-l<n 



Ji...i„_i 



det {YX*) (Ti • • • cr„_i det ( ( Y* A (0) X 

n-1 



l...(n-l) 



det {YX*) cTjCaCi (uq, ^' (0) wq) 7^ 0. 
7=1 ^ ^ 



This completes the proof. □ 
We now prove Theorem 13.1 



4.2. Proof of Theorem 13.11 Proof. We begin by noting that our hypotheses 
imply that statements (ii), (iii), and (iv) of Theorem 12.11 are true. In particular, 
statement (iii) implies that there is exactly one convergent Puiseux series for the 
Ao-group given by 

A„ (e) = Ao + AiC^" + (^"^") ' ^ = 0' - 1' 

k=2 
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where C = and Ai ^ 0. Then by weU known results [H §6.1.7], pi §11.1.2] there 
is a convergent Puiseux series for the corresponding eigenvectors of the form 

Xh (e) = xa + '^Xk (^C'^e™) , ^ = 0, m - 1, 

where C = e^' and xq is an eigenvector oi Aq ^ A (0) corresponding to the eigenvalue 
Aq. Now if we examine the proof of (ii)=>(iii) in Theorem 12.11 we see by equation 
(|4.3p that A™ = £~^, where £,„ is given in equation (|4.2p in the proof of (i)=>(iii) for 
Theorem 12. II Thus we can conclude that 

,™ ^ f det(A/-A(£))|(,.,)^(o,,„) 



Choose any mth root of {vm, AiUi) and denote it by (vm, AiUi)^^"^ . Applying 
this mth root to both sides of the above representation for A™ we find that 



/ \ 1/™ 

^det(A/-^(e))|(,,A)=(o,A„) 



v 



^ dct(A/-A(E))|(^,A) = (O.Ao) 



for some r depending on Ai with < r < m — 1. But then 



X^) = Xh+r (e) = Ao + ^ Afc (c'^+'e™ ) 

fe=i 

1 °° fe 
= Ao + (AiC'')(c'e-) ,/i = 0,...,m-l 

fe=i 

is still the same collection of eigenvalues for A (e), since = ^[{h+r) mod m] ^ 
still a convergent Puiseux series with corresponding eigenvectors 

oo 

Xh (e) = Xh+r (e) = Xq + ^ u'^^'^e" 



fe=i 



Hence without loss of generality we may assume that, for the Puiseux series repre- 
senting the eigenvalues of A (e) given in p.7p . 



Ai 



l/m 

-§^det{XI-A{e))\(,^x)=io,Xo) ' 

' ^ dct(A/-A(e))|(^,A)^(o,AQ) 



Next, we wish to prove that we can choose the perturbed eigenvectors (|3.8p to 
satisfy the normalization conditions (|3.10p . For us to do this we must first show that 
{vi,Xo) 7^ 0. This is done by noting that Theorem l2.1l fiv) tells us the Jordan normal 
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form of A (0) — Aq corresponding to the eigenvalue Aq consists of a single m x m 
Jordan block. Hence the kernel of Aq — XqI is one-dimensional so that there exists a 
Ci ^ such that xq = ciUi, since xq and ui are eigenvectors of Aq corresponding to 
the eigenvalue Ao( the latter following from Proposition lA.il (ii)). Then by this and 
Proposition lA.il (i) we have {vi,xo) — ci ^ 0. 

Now since in a neighborhood of the origin Xh (e) is an eigenvector of A (e) corre- 
sponding to the eigenvalue Xh (e) and Xh (s) is continuous at the origin with Xh (0) = 
Xq then in possibly a smaller neighborhood of the origin [vi, Xh (s)) ^ and x^ (e) = 
(vi'xhie)) eigenvector of A{e) corresponding to the eigenvalue Xhis)- But 

then xo{z"^) = ^^^^^^J|^2.^^ is analytic at the origin since 2:0(2:™) and {vi,xq {z"^)) 
are. Hence we have the series representation xo{z"^) = xq + J2'k'=i^^'^ ^ which is 
convergent in a neighborhood of the origin. This implies x^ (e) = x^ (^(^e^^"^^ = 

^"'"Efe^i ^ (^C^s^^"^^ for h — 0, m— 1. But xq (z™) is analytic at the origin and 

hence so is {vi,x^{z"')) with 1 = (vi, (;;^|^iJrj)) = (z™)) EfeLi {vi,xj:)z''. 

This implies (ui,5Jo) = 1 and {vi,x^.) = 0, for fc > 1 and {vi,xii{e)) — 1 for 
h = 0,...,m— 1. Combine this with the fact that x^ (e) is an eigenvector of A{e) 
corresponding to the eigenvalue A;i (e) and we have proved the statement that we 
choose the perturbed eigenvectors (|3.8p to satisfy the normalization conditions p.lOp . 
Note that we have also proved that for those perturbed eigenvectors that satisfy the 
normalization conditions we have {vi,xq) = 1 and (vi,Xk) = 0, for fc > 1. 

We are now in a position to show that Ai, A2, ... and xq,xi^ ... are given by the 
recursive formulas (|3.1ip - (|3.13p . First some preliminaries. We begin by using the 
equality 

= (e) - Ao (e) /) xq (e) 

/ 00 \ 00 

) k=Q 

- Afe/ ) Xs-k 




in a neighborhood of the origin, where Ak_ = if ^ N. This equality implies that 
Efc=o i^-!^ ~ Afe/j Xs-k ~ 0, for s > 0, and hence 



(Aq - XqI) Xs = - J2 (^^ " ^kA Xg-k, for S > 1. 
k=l ^ " ^ 

Next, denote 5* by the space spanned by the vectors {Uei\ 2 < i < n}. We will 
now show that Xg € S for s > 1. First notice that since U is invertible and is an n x n 
matrix then its columns, Uci for 1 < i < n is a basis for C"^^. Hence we have for 
each s with s > 1, Xs — E"=i (^iUei for some complex scalars Next, we note that 
by the facts {vi, xo) = 1 and {vi,Xk) — 0, for fc > 1 as we just proved, the definition 
vi = {U^^) ei, and (wi, Uci) — Su for 1 < i < n, this implies that — {vi,Xs) — ci 
for s > 1. Hence Xs & S for s > 1. 

Now from Proposition lA.ll (iii) we know that A {Aq — XqI) acts as the identity 
on S* = span{Uei\ 2 < i < n}. Hence since Xg (z S for s > 1 we have Xg = 

A {Aq - Aq/) Xs, for S > 1. 
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Finally, for s ~ we have Xq = ui. This follows from that facts that we have 
already proven that there exists a Ci 7^ such that xq = CiUi, the normalization 
condition (wi, xq) = 1, and (wi, Uei) = 1. Thus we have proved that 

xq = Ml, Xs ^ A {Ao - XqI) Xs, for s > 1, (4.4) 

which is the first key step in this proof. 

The next key step in this proof is the following lemma. 
Lemma 4.1. For all s > the following identity holds 



{Aa - Ao/) Xs = < 



E 

1=0 



/ PsjUi, for < s < m 



ni s~m j \_ m \ 

Y,Ps,iU.i- J2 J2 J2 Pj,kA''AiXs-j-ir,i, for s > m 

i=0 j=0 k=0 1=1 



(4.5) 



where we define uq — 0. 

Proof. The proof is by induction on s. The statement is true for s = since 
{Ao — XqI) xo = = pOjOWq. Now suppose it was true for all r with < r < s for 
some nonnegative integer s. We will show the statement is true for s + 1 as well. 

Suppose s + 1 < m — 1 then {Aq — XqI) Xr — X]I=o^''"^«"^ tor < r < s and 
we must show that {Aq — XqI)xs+i = J2i=o Ps+i,i'^i- Well, first for 1 < r < s 
we have a;^ = A {Aq — XqI) Xr = J^l^o Pr,iAui = Pr,iAui by the first key step 

(|4.4p . Second, by Proposition lA.il (iv) we have Am = Ui+i for 1 < i < r and hence 
Xr = J2i=iPr,i'^i+i- Hcuce wc have 



s+1 / X 

(^0 - Ao/) xs+i = - E [Aj^-Xkl)x 
fc=i ^ " ^ 



s+l-fc 



s s+1 — 

Xs+iUi + J2 E ^kPs+l-k,iUi+i 
k=l 1=1 

s+1 s+l-fe 

J2 J2 ^kPs+l~k,iUi+l, 
k=l i=0 



the latter equality holding since Ps+i-fc.o = for s + 1 — A; > and po,o = 1- Now 
we rearrange the sum using the identity Efc=i Ei=o = Ei=o Efe=i ' ^k,i and 

use the fact that J2kt\^' ^kPs+i-k,i = Y}k=i^~'^'^^^^^ ^kP{s+i)-k,[i+i)-i = Ps+i.i+i 
for s > i > by Proposition lB.il (vi), to get 

s /s+l-i \ 

(^0 - Ao/)a;>;+i = E E ^kPs+i-k,i ) "i+i 



i=0 \ k=l 
s+1 

= E Ps+l,iUi. 

1=0 

Therefore the statement is true if s + 1 < m — 1. 

Now suppose that s + 1 > m. The proof is similar to what we just proved. First 
we get for 1 < r < s, = A [Aq — XqI) Xg, for s > 1 by the first key step (|4.4p . 
Aui = Ui+i for 1 < i < m — 1 and Aum = by Proposition lA.ll (iv), and Pr^ = 
by definition since r > 1. Hence these together with the induction hypothesis for 
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1 < r < s that (14. 5p is true implies that 

r 

J2 Pr,iU'i+i, for < r < TO - 1 

1=0 

m— 1 r — 'm j \_~m~\ 

J2 Pr,iUi+l - J2 J2 J2 P],k^''^^AlXr-j-lm, for r > TO 
i=0 j=0 k=0 1=1 

Thus from this representation of Xr for 1 < r < s we get 

{Ao - \oI) Xs+i = - [Ajl - Afe/ Xs+i-fc 
fe=i ^ 

L^J .+1 

= — X] AiXs+l-lm + J2 ^kXs+l-k 
1=1 k=l 

— ^ J2 AiXs+i-lm + J2 ^kXg+l-k + ^kXg+l-k 

1=1 s>s+l-k>m m>s+l-k>Q 

= - I] AAiXs-lm + Ps+l-k,iUi+l- 
1 = 1 s>s+l-k>rn \ i=0 

s+l-k-m j ["'^''m'' '\ 

Y Y Y Pj.jA'+MzXs+l-fc-j-im 

j=a 1=0 1=1 

/s+l-k 

+ Y Y Ps+l"fc,iUj+l 

m>s+l-fc>0 V i=0 

I £±i I 

L m J s+1— mm— 1 

= - Av4ia;s_im + X] S ^kPs+i~k,iUi+i 

1 = 1 k=l i=0 

I B+l-fc-3 I 

s+1— ms+1 — fc — m 3 I m j 

- Y Y Y Y XkPj,i^'''^^ AiXs+i-k-j-im 

k=l j=Q i=0 1=1 

s+1 s+l-k 
+ Y Y XkPs+l-k,iUi+l. 

k>s+l — m i=0 

We use the sum identity Yltl+i-m Yt=o''' = YT=o^ Yltl7i-m «fe,' and the 

identity YltV ^kPs+i-kA = Yk=i^~^''^^^^^ \kP(s+i)-k,(i+i)-i = Ps+i,»+i for « > 
I > (which foUows from Proposition [RT] (vi)) to get 

I s+1 1 

L Tlx J m— ls+1 — m 

{Ao - XoI)Xs+l^ ~ Y AAiXs^lm+ Y Y XkPs+l~k,iUi+l 
1=1 i=0 k=l 

I s+l-k-j I 

s+1 — m s+l — k — 'm J L m j 

- Y Y Y Y AfePj,iA'+M;Xs+i_/c-j-;m 

k=l j=a i=Q 1 = 1 

m — 1 s+1— 2 

+ Y Y XkPs+l-k,iUi+l 

i=0 /e>s+l — m 

L^^J m-1 /s+l-i \ 

= - Y AAiXs-lm + Y { Y XkPs+l-k,i I Wi+1 
1=1 i=0 \ k=l J 
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I s+l-fc-j I 

k=l j=0 1=0 1 = 1 

= - ^AXs-lm + J2 Ps+l.i+lUi+1 

1=1 i=0 

I I 

s-\-l~ni s-\-l~k — m J \_ m j 

- J2 J2 J2 \kPj,i^''^^AlXs+l-k-j-lm 
k=l j=0 1=0 1 = 1 

= - J2 ^AiXs~lm + J2 Ps+l,iUi 

1=1 i=a 

I s+l-fc-j I 

s+1— ms+1 — fc— m J \_ rn J 

- J2 Y J2 J2 XkPj,2^'^^AiXs+l-k-j-lrn- 
k=l j=0 1=0 1 = 1 

Now though we use the sum identities 

s+l — m s+l — fc — m j s+1— m— 1 J — m j 

E E E ofc.i,* = E EE ofcj.i 

fc=l j=0 i=0 j=0 k=l 2=0 

s+1 — m— 1 i s + 1— j — m 

= E E E 

j^O i^O A.— 1 

s-^l—m—1 — m — 1 s+1— J — m 

= E E E 

i—O j—i k—l 

s+1— m— 1 s+1 — m 

= E E E «fcj,i 

i=a q=i+l j+k=q 

i<j<s+l — m — 1 
l<A;<s+l— J— m 

s+1— m— 1 s+1 — m q — i 

= E EE 0'k,q-k,i 

1=0 q=i+l k=l 

and let akj.i = E XkPj^iJ^''^^ AiXs+i-k-j-imi with the identity ^kPq-k.i = 

1=1 k=l 

g-(j+l) + l 

E XkPq-k,(i+i)-i = Pq,i+i for (7 > i + 1 > 0(wliich follows from Proposition lB.il 

k=l 

(vi)) to conclude that 

s+1— m s+1 — fc— m j 

E EE o.k,j,i 

k=l j=0 2=0 

s+1— m— 1 s+1 — m q— 4 

E EE 0.k,q-k,i 

1=0 q=i+l k=l 







m 


E ^Aixs- 


-Im 4 


- E Ps+l,iMj 


1=1 




i=0 


L^J 




m 


E ^AiXs- 


-Im ~1 


- E Ps+l,iWi 


1=1 




i=0 






m 


E ^AiXs- 


-Im -1 


- E Ps+l,iMi 


1=1 




i=0 


s+1— m— 1 s+1- 


-m ij- 





L 2^ 2^ AkPq-k,iI\ AlXs ~\-l — q—lm 

^0 g^i+1 fc^l 1^1 
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■ J2 ^AXs-lm + J2 Ps+l,iUi 
1 = 1 1=0 

I I 

s+1— -m— 1 s+1— m L m J / q—i 

E E E E AfeP,-fe. ) A^+M,x,+i„,_i„ 

i=0 i3=i+l i = l \fc=l 



1=1 



i=0 



s+l-m-1 s+l-m ["^m 'J 

- E EE PqA+l^-'^^AlXs+l-q-lm- 

1=0 q=i+l 1=1 

Next we use the sum identities 

s+1— m— Is+l — m s+1 — m<?— 1 s+l—mj — 1 s+1 — m j 

E E s.i+i = E E = E E «j.fc+i = E E »i,fc 

1=0 q=i+l q=l i=0 3 = 1 fc=0 j=l k=l 

together with the fact that o = for j > and po o = 1 to conclude by letting 

= E Pq.i+l-l^''^^ AlXs+l-q-lm. that 
1=1 



[Ao - Aq/) x^+i 



1 = 1 i=0 

s+l-m j L"^™ 

- E E E Pj.kJ^^AiXs+i-j-im 

j=l k=l 1=1 

I s+l-j I 

m s+l-m J L TTi J 

= J2 Ps+l,iUi - E E E Pj,kA-''AlXs+l-j-lm- 
i=0 j=0 k=0 1=1 

But this is the statement we needed to prove for s + I > m. Therefore by induction 
the statement (|4.5I) is true for all s > as desired. □ 

The lemma above is the key to prove the recursive formulas for As and Xs as 
given by (|3.1ip - (|3.13p . First we prove that Xs is given by (|3.13p . For s = we have 
already shown xq — ui = po,oUi. So suppose s > 1. Then we use the above lemma 
identity (|4.5p together with the identities given in (14. 4p . the fact that Aui = Ui+i for 
1 < z < TO — 1 and — (see Proposition I A. ll (iv)), and Ps.o = by definition 
since s > 1 to conclude that 



Xs = A {Ao - XqI) Xs 



= < 



i=0 



m—1 



"y^^Ps^iAuj, for < s < TO — 1 



1=0 

s-m j L~sr"J 



Y^PsA^Ui - E E E Pj,fcA''+lAiXs-j-im, for S > TO 



j=0 /c=0 1=1 
s 

^Ps,jUi+i, ifO<s<m-l 



1=0 

s-m j L~sr"J 



Ps,iUi+i- E E E P].k^^'^^AiXs-j-im,'^is>m 

4=0 j=0 fc=0 1=1 
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This proves that is given by (|3.13p . 

Next by using the above lemma identity (|4.5p again we will prove that is given 
by p. lip and (|3.12p . We start with s = 1 and prove Ai is given by p. lip . To being 
we note that by Proposition lA.ll (ii) we have {Aq — XoI)*Vm = 0. Then using the 
above lemma identity (|4.5p for s — m, the biorthogonality relation (vi, Uj) = dij from 
Proposition lA.il (i), and the fact pm.m = A™ from Proposition lB.il (iv), we get 

0= {{Aq - XoI)*U,n,Xrn) 
= (w„i, {Ao - XoI)Xm) 
( 



\ 



ra ra—m j l_ m J 

i=Q j=0 k=0 1=1 



= (I'm, A"Um - AiUi) 
= A™ - [Vm,AiUi) , 

SO that A™ — (vm, AiUi). Now in the beginning of our proof we chose a mth root 

of {vm,AiUi) and denoted it by (wm, ^iwi)^^™- We showed that for this same mth 
root we had 



/ 



l/m 



Ai 



^det(A/-A(£))|(,^;,)^(0,;,„) 

dct(AJ-A(e))|(e,A) = (o,Ao) 
m! 



This means that 



^det (A/- A(£))|(,,A)=(o,Ao) 



/ a" 
( Ml 



■ dct(A/-A(£))|(e,A) = (0,Ao) 



A™ = {vm,Aiui) 



and by taking that mth root of {v„n AiUi) on that chain of equalities we conclude 
that formula (|3.1ip is true. 

Finally, suppose that s > 2. Then using the fact (Aq — Ao/)*Wm = by Proposi- 
tion |X]T](ii), the above lemma identity (|4.5p for m + s— 1, the biorthogonality relation 
(vi^Uj) — Sij from Proposition lA.l] (i). the identities A*vi = 0, and A*Vi = Ui-i, 
for 2 < i < m from Proposition lA.ll (v), the sum identity J2'j=o^k=o^j-k = 



J2k=o J2''i= k and the relation rm+s-i,m = Pm+s-i,m - wA™ ^A^ from Propo- 
sition IB. II (ii) , we get 



= ((Ao - \oI)*Vra,Xra+s-l) 
= {Vm, (^0 - XaI)Xrn+s-l) 



I 



(m+s — 1) — m J L 



(m+s — 1)— J — /r) 



^m 7 / ^P7n+s — l,i'^i 
i=0 

s-1 j 



3=0 



k=0 1=1 



Pm+s — l.m 

ZZ^'J''^ i^*)''^ra, AiXm+s-l-j- 
j=0 fc=0 \ i=l 
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s-1 s-1 
■i=0 j=i 
m — 1 s — 1 

2 = J=I 



"m-i, E ^fcX„i+s-i-i-fem I , for 2 < s < m 
fc=l 



^fc^m+s — 1— J — /cm 

fe=l 



, for s > m 



And therefore using this equahty and the fact Ai 7^ we can solve for and we will 
find that it is given by (|3.12p as desired. This completes the proof. □ 
Finally, we can now prove the corollaries. 

4.3. Proof of Corollaries 13.21 and 13.31 Proof. Both corollaries follow almost 
trivially now. For if we examine the recursive formulas in Theorem 13. II p.lip - p.l3p 
we see that Afc and Xk require just Aq, . . . , Aj^ m+fc-i j . Furthermore, to get up to the 
second order terms we just explicitly calculate A2 and xi, X2- The only nontrivial 
aspect is the calculation of A2. But this follows from the series representation of 
/ ^) (|4-f P ^"^d Ao (e) p.7p which, for e in a neighborhood of the origin, satisfies 



= / (e, Ao (e)) a^e + aom (Ao (e) - Aq)" + 



E 



a,ye* (Ao (e) - Aq)^ 



(ij)^{(Oj):j<m} 



k—m 

= (Oio + aOmPm.m) £ 



E E 



i+j=2, iJGN k=j 



for m = 1 



{aOmPm+l,m + aOm+lPm+l,m+l + OllPl,! + a20POo) £ 

{aOmPm+l.m + aOm+lPm+l,m+l + ailPl,l) , for TO > 1 



-O 



= (flio + aomA™) e 

_^ I (aom'TiA"^^ A2 + flom+i A™'^^ + an Ai + 020) , for m = 1 
I (aommA™"^A2 + ao„i+iA™+^ + anAi) , for m > 1 

(see the proof of (i)^(ii) in Theorem (12. ip and Proposition lB.il (iii)-fv)). This allows 
us to solve for A2 yielding 



A2 = 



for TO > 1 



-(aom+l A'"+^+aii Ai+a2o) 
aommA"^ 

^^^^(O.Ao)+Ai^(0,Ao) + l|^(0,Ao)) 



'1 (m+l)! aA-i 



-|,Ai („+!)! 3Arr. + l 



|^(0,Ao))^ 



for m — 1 



:(0.>o)+Ai^(0,Ao)) 
™^r" (iT|;;*4(o,AQ)) 

And therefore the corollaries are true. □ 

Appendix A. 
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Let Aq and U he n x n matrices with the latter being invertible such that 



Jm (Ao) 






Wo 



(A.1) 



where Wq is a (n — m) x (n — m) matrix such that Aq is not one of its eigenvalues. 

Define the vectors . . , wi,. . . , Vm S C"^"^ as the first m columns of the 
matrix U and (C^~^) , respectively, i.e., 



And define the matrix A € 
k:=U 



Ui : = Uci, 1 < i 


< TO, 


Vi : = ([/- 




1 < i < m. 


-nxn 






r Jm (0)* 








(Wo- 





(A.2) 



where is the (n — m) x (n — m) identity matrix. 

Proposition A.l. The following properties hold: 
(i) The two sets of vectors {ui : I < i < m} and {vj : 1 < i < m} are 
biorthogonal, i.e., 



{Vi,U 



where 6ij is the Kronecker delta. 

(ii) The vector Ui is an eigenvector of A^ corresponding to the eigenvalue Aq 
and Vm is an eigenvector of Aq corresponding to the eigenvalue Aq. 

(iii) We have A {Aq - Aq/) Vex = 0, A (Ao - Aq/) Vci = Uci, for2<i<n, 

(iv) Aum = 0, Aui = Wj+i, for 1 < i < m — 1, 

(v) A*vi = 0, and A*Vi = Wj-i, for 2 < i <m. 

Proof, i. The biorthogonality follows from fact {ei,ej) = Sij since 

(^{U^^y ei,Uej^ = {ei,ej) . 
ii. This follows from the definitions Ui = Uei, Vm = {U~^Y e^ and the 



facts 



Jm (Ao) 






Wo \ 



ei = AqCi, 



Jm (Ao) 






Wo _ 



m. This follows the fact that Jm (0)* Jm (0) = dmp[0, /m-i] since then 
A(Ao- Ao/)t/ = AC/ 



■ Jm (0) 






(Wo - Ao/n-m) . 



r (0)* 






\ Jm (0) 






(Wo - Xoln-m) ^ . 




(Wo - Xoln-m) 



= u 
= u 



iv. From the definition Ui = Uei, 1 < i < m it suffices to prove that 
U~^AUem ~ and U^^AUci = Cj+i, for 1 < i < m — 1. But this just follows from 
the definition of A and Jm (0)*. 
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V. From the definition Vi — {U ^)* e^, 1 < i < m it suffices to prove ttiat 
U*h* (C/-i)*ei = and U*h* (C/"i)*e, = e,_i, for 2 < i < to. But this foilows 



from the fact 



' Jm (0) 








(Wo - Ao/„-m) ^ 





Remark A. 2. There is another way to define vi, ■ ■ ■ , Vm o,nd A without needing 
to put Aq into the block representation iA.l\) . Start with any Jordan basis ui, . . . , 
Um with generator Um for Aq corresponding to Xq, i.e., {Aq — XqIY Um = Um-i, for 
i — 0, . . . , TO — 1. Then there exists unique vectors vi, ■ ■ ■ , Vm G C"^^ such that for 
Uq := • • • \um] and Vq := • • • \vm] we have 



Pa 



where P\g € ^nxn eigenprojection of Aq corresponding to Aq, i.e., the projection 

whose range is the algebraic eigenspace of Aq corresponding to Xq. Define the matrix 



A UqJ^ (0)* Vq* ~ Sq, 

where Sq G ig the reduced resolvent of Aq at Xq '1, ^2.6.3 Reduced Resolvants, 

pp. H-15, Eq. (6.6)]. The matrix Sq has the following integral representation 

where R (A) = [XI — Aq) ^ is the resolvent of Aq and the integral is over any simply 
closed positively oriented piecewise smooth path on which the resolvent is defined and 
which its interior contains no other eigenvalue of Aq except Xq. If we extend Ui, 
. . . , Um to a basis ui, . . . , Um, Um+i, ■ ■ ■ , Un for C"^^ and define the matrix U := 
• • • |Mm|Mm+i| ' ' ' \un\ then Aq would havc the block representation \A.1\) and we 
would find that Vi = {U~^^ a and A is equal to the right hand side of \A.^) . 

Appendix B. 

Foilowing 7, §9.32 pg. 414] we define the poiynomials pj^i — pj^i (Ai, . . . , Xj-i+i) 
in Ai,. . . , Aj_i+i, for j > * > 0, as the expressions 

Pos) ■ = 1, Pjfi ■= 0, for j > 0, 

Pj,, (Ai, . . . , Aj_,+i) : = 22 n Asg, for j > i > 0, 

l<Se<J-«+l 

and the polynomials rk+i,k = rk+i,k{Xi, . . . , A;) in Ai,. . . , A;, for A:,; G N, as 

k 

rk+i.k{Xi, Xi) := 0, rk+Lk{>^i, . . . , A/) := 2_, 11 ^s^, for / > 1. 

si + --- + Sk=k+l 



Assume that X^jli '^j-^'' ^ convergent Taylor series. 

Proposition B.l. The polynomials defined above have the following properties: 
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(i) J:Pj,iZ' = E A,^^ ] ,forj>i>0, 



have 



(ii) rk+Lk = Vk+Lk - fcAi"^A;+i, for k,l&n, 

(iii) = \j, for j > 1, 

(iv) Pj,3 = /or j > 0, 

(v) pj+ij = jX{~^X2, forj > 0, and 

j-i+l j-i+l 

(vi) J2 XqPj-q,i-l = E Pq,lPj-q,i-l "^Pj.t, forj > 1 > 0. 
q=l q=l 

Proof, i. For i = we have ^Ejli Aj-^'^^ = 1 = Ejli Pjfi^^ ^ ^'^'i for i > we 



OO OO 



EvM =E-E HA, 



yS\-\ \-Si 



Sl=l 5i=l 

I 



E 



E HA., 



siH |-Si=j 

\l<Se<j-»+l 



M. Let fc, / e N. Then by definition of Pk+i,k we have 

fc 



E UK 



si + ---+Sk=k+l s^-- ' 

l<Sg<l+l 

E 

siH hSfc='i:+i 

l<Sg<J+l 

3g^{l,...,k} such that s^— /+1 
k 



n A.,+ 

o=l 



E 



.s-i+---+sfc=fe+; 

l<S(j<l+l 

such that Sg—l+l 



k 

n A. 



^ArAm+ E nA« 



vfe-l 



siH |-Sfc=fe+i 

l<s„<i 



= fcAi Ai+i+rfe+i,fe. 
m. By the definition of pj i we have Pji = ^ si=j no=i As = Aj , as 

desired. 

it;. By the definition of pj,j we liave po,o = 1 for i = and for j > we have 
Pj,3 = J2si+-+sj=j Yll=i Asg = 11^=1 Ai = Aj. 

l<Se<l 

u. By the definition of Pj+ij for j > we have 



Pj+K 



E n A., = EAr'A2=jAr'A 



siH hSj=j + l 

l<s„<2 
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vi. Let j >i > 0. Then by part (i) we have 



E 



I \ 



= E 



ji>i-l 

i2>i 



oo /j-i+l \ 

= E ( E Pj-l'i-^Pl^ ) 
j=i \ g=l / 



implying together with part iii of this proposition that pj^i = '^tJi^^ Pj-q,i-iPq,i 



Eg=i Pj-q,i-i\ for j > t > 0. □ 
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